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ABSTRACT 

Pulsed high energy radiation from pulsars is not yet completely understood. In this paper, 

■ we use the 3D self-consistent annular gap model to study light curves for both young and 

millisecond pulsars observed by the Fermi Gamma-ray Space Telescope. The annular gap 

O ' can generate high energy emission for short-period pulsars. The annular gap regions are so 

large that they have enough electric potential drop to accelerate charged particles to produce 
gamma-ray photons. For young pulsars, the emission region is from the neutron star surface 
to about half of the light cylinder radius, and the peak emissivity is in the vicinity of the null 
charge surface. The emission region for the millisecond pulsars is located much lower than 
that of the young pulsars. The higher energy 7-ray emission comes from higher altitudes in the 

■ magnetosphere. We present the simulated light curves for three young pulsars (the Crab, the 

Vela, the Geminga) and three millisecond pulsars (PSR J0030+0451, PSR J02 18+4232, PSR 

CN| ■ J0437-37 15) using the annular gap model. Our simulations can reproduce the main properties 

C") I of observed light curves. 

■ Key words: Radiation mechanisms: non-thermal - Pulsars: individual: Crab, Vela, Geminga, 
O ! J0030+0451, J0218+4232, J0437-4715 - Gamma-rays: stars. 



1 INTRODUCTION 

Pulsars are fascinating astronomical objects in the universe. After more than 40 years since the discovery of the first pulsars, their pulsed 
non-thermal emission has not been completely understood due to the insufficient knowledge about the global acceleration electric field and 
particle dynamics in the magnetosphere. 

High energy emission (e.g. 7-ray emission) from pulsars takes away a significant fraction of the rotational energy. Several space tele- 
scopes were used to observe high energy emission from pulsars. Six 7-ray pulsars and three candidates were discovered by the Energetic 
Gamma Ray Experiment Telescope (EGRET) (Thompson et al. 1999). One young pulsar (PSR B 1509-58) was detected only up to 10 MeV 
by the Imaging Compton Telescope (COMPTEL) (Thompson 2001). Pulsars detected by EGRET show double peaks with bridge emission. 
Thanks to the launching of Astro-rivelatore Gamma a Immagini LEggero (AGILE) and Fermi Gamma-ray Space Telescope (FGST), more 
than forty new 7-ray pulsars have been discovered in the last year, including gamma-ray only pulsars and a new population of millisecond 
pulsars (Pellizzoni et al. 2009, Abdo et al. 2009b, Abdo et al. 2009c). With these data, we have good opportunities to study the open ques- 
tions, e.g., the location of the emission zones, possible particle dynamics in the magnetosphere and so on. Theory on non-thermal high energy 
emission from pulsars should be significantly improved in the coming years. 

Physical and geometric magnetosphere models have been proposed to explain pulsar's 7-ray radiation, as we will summarize below, 
differing on the acceleration region of the primary particles and the mechanism for the production of the high energy photons. 

The first is the polar cap model (Daugherty & Harding 1982, Daugherty & Harding 1994, Daugherty & Harding 1996). The acceleration 
region is located in the vicinity of neutron star surface up to tens of kilometers near the magnetic pole. The observed phase-resolved spectra 
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of 7-ray pulsars can also be modeled. This model favors a small inclination angle, e.g., a nearly aligned rotator. The 7-ray emission region is 
so close to the neutron star surface that high energy 7-ray photons could be absorbed in the strong magnetic field (Lee et al. 2010). 

The polar cap model was combined with the slot gap model recently to explain the pulsar 7-ray radiation (Muslimov & Harding 2003, 
Muslimov & Harding 2004). A 3D model of optical-to-7-ray emission from the slot gap was developed to study the Crab pulsar's light curve, 
phase-averaged and phase-resolved spectrum (Harding et al. 2008). Hirotani (2008) demonstrated that the slot gap model reproduces at most 
20% of the observed fluxes because of the small trans-field thickness. The two-pole caustic model is proposed by Dyks & Rudak (2003). The 
gap is thin, confined to the surface of the last open field lines, and extends from both polar caps to the light cylinder. Double peak light curve 
with a large peak separation (e.g.. Vela) can be well reproduced by this model. But this model has difficulties in explaining the light curves 
with small peak separations (e.g., B 1706-44 and B 1055-52). Frackowiak & Rudak (2005) presented spectral and light curve properties of 
gamma radiation obtained by numerical modeling of some pulsars, i.e. PSR J02 18+4232, PSR J0437-4715 and PSR B 1821-24. 

The outer gap model is another excellent model to interpret 7-ray emission from pulsars (Cheng, Ho & Ruderman 1986a, Cheng, Ho & 
Ruderman 1986b, Chiang & Romani 1992, Romani & Yadigaroglu 1995, Romani 1996, Zhang & Cheng 1997, Cheng, Ruderman & Zhang 
2000, Zhang et al. 2004, Lin & Zhang 2009). The classical outer gap starts at the null charge surface (inner boundary), ends at the light 
cylinder (outer boundary); it is further bounded by the last open field line (lower boundary) and a layer of electric current (upper boundary). 
However, Hirotani et al. (2003) argued that the position of the inner boundary could be shifted towards the neutron star surface because a 
current at nearly the Goldreich- Julian rate (Goldreich & Julian 1969) is injected to the outer gap. Tang et al. (2008) used this modified outer 
gap to study the multifrequency phase-resolved spectra of the Crab pulsar. 

Watters et al. (2009) have simulated a population of young pulsars and computed the beaming pattern and light curves for the three 
models: the polar cap model, the two-pole caustic (slot gap) model and the outer gap model. Venter et al. (2009) presented light curves of 
millisecond pulsars from 3D emission modeling, in the geometric context of the polar cap, the outer gap, and the two-pole caustic models. 
They found that most of the light curves are best fit by the two-pole caustic and the outer gap models, which indicates the presence of narrow 
accelerating gaps limited by robust pair production - even in these pulsars with very low spin-down luminosities. 

Force-free relativistic MHD were used to study the time-dependent evolution and dynamics of pulsar magnetospheres for either aligned 
or oblique magnetic geometries (Contopoulos et al. 1999, Komissarov 2006, Spitkovsky 2006, McKinnqey 2006, Timokhin 2006 and Gruzi- 
nov 2007). Recent efforts on theoretical understanding the high energy emission from pulsars can be found from Bai & Spitkovsky (2009a; 
2009b). They modeled gamma-ray pulsar light curves for the two-pole caustic model (slot gap model), the outer gap model and the "separa- 
trix layer" model using the more realistic magnetic field taken from 3D force-free magnetospheric simulations. Their separatrix layer model 
might be associated with the current sheet beyong the light cylinder. Their simulated results indicate that the separatrix layer model can best 
reproduce the observations. 

The annular gap model we are developing is originally proposed by Qiao et al. (2004) and Qiao et al. (2007). The gap region is located 
between the critical field lines0 and last open field lines, and extends from the neutron star surface to the light cylinder. The region for high 
energy emission in the annular gap model is concentrated in the vicinity of the null charge surface, i.e., an intermediate emission height, 
different from the outer gap model. The annular gap has a sufficient thickness of trans-field lines and a wide altitude range for particle 
acceleration. The role of the annular gap depends on the mono-polar voltage in the annular region, so it is favorable for short period pulsars. 
This model combines the advantages of the polar gap, the slot gap and the outer gap models, and works well for pulsars with short spin 
periods. It is a promising model to explain high energy emission from young and millisecond pulsars. 

In this paper, we focus on the 7-ray light curves for young and millisecond pulsars in the framework of the annular gap model. In § 2, 
detailed physics and geometry of the annular gap are introduced. The results of simulated photon sky-maps and light curves for young and 
millisecond pulsars are presented in § 3. Discussions and conclusions are given in § 4. 

2 THE ANNULAR GAP MODEL: PHYSICS AND GEOMETRY 

The open field line region of a pulsar magnetosphere can be divided into two parts by the critical field lines. One is the annular region, which 
is between critical field lines and last open field lines. The other is the core region, which is within the critical field lines. Taking an anti- 
parallel rotator as an example, the radii of the core polar region (r- CO rc) and the full polar cap region (r p ) are r cor c = (2/3) 3//4 i?(f2i?/c) 1//2 
and r p = i?(f2i?/c) 1//2 , respectively (Ruderman & Sutherland 1975), here R is the pulsar's radius, D. is the angular spin frequency. The 
radius of the annular polar region is r an n = r p — r COIC = 0.26i?(f2i?/c) 1 ^ 2 . If pulsar's spin period is smaller, the annular radius (r ann ) is 
larger. Therefore, the annular acceleration region is more important for those pulsars with a short spin period (e.g., millisecond and young 
pulsars) and is negligible for older pulsars with a large spin period. 

If pulsars are bare quark stars (Xu 2002, 2005), then two kinds of acceleration regions could be formed, namely, the inner vacuum core 
gap and inner vacuum annular gap. If pulsars are neutron stars with enough surface binding energy, the inner vacuum core gap (Ruderman 
& Sutherland 1975) can be formed for an inclined rotator with a magnetic inclination angle a > 90°; while the inner vacuum annular gap 
(Qiao et al. 2004) can be formed only for an inclined rotator with a magnetic inclination angle a < 90° . If the binding energy of a neutron 
star's surface is not high enough, both negative and positive charges will flow out freely from the surface. Then two kinds of acceleration 

1 The critical field line is across the intersection of the null charge surface and light cylinder. 
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Figure 1. The coordinate systems for an oblique rotator. In the laboratory frame O — XY Z, the Z axis is aligned with the rotational axis f2 

of a pulsar. The magnetic frame O — X'YZ' is generated by rotating the O — XY Z coordinates around Y axis by an inclination angle of a, the Z' axis is 
aligned with the dipolar magnetic moment /i. Point "s" is an emitting source on an arbitrary magnetic field line. The r s is the altitude of the radiation source 
from the neutron star center. The polar angle and azimuthal angle in laboratory polar coordinates are denoted as ©,</>, while the polar angle and azimuthal 
angle in the magnetic polar coordinates are 0, tp. 

regions (annular and core) could be formed. Particle acceleration is more effective in the annular gap region. The annular acceleration region 
extends from the pulsar surface to the null charge surface or even beyond it. This leads to a fan-beam 7-ray emission, which is suitable to 
interpret the observed light curves and the broad-band emission (Qiao et al. 2007). 



2.1 Coordinates for the Annular Gap 

As shown in Fig.l, in the coordinate O — XY Z of the laboratory frame, the Z axis is aligned with the rotational axis Q of the pulsar. The 
magnetic frame O — X'YZ' is generated by rotating the O — XY Z coordinate around Y axis by the inclination angle a, the Z' axis is 
aligned with the dipolar magnetic moment /i. The two vectors and /i locates in the plane O — XZ, which is called the fl — /x plane. The 
polar coordinate associated with O — XY Z and O — X'Y Z' are called the laboratory polar coordinate and the magnetic polar coordinate, 
respectively. The polar angle and azimuthal angle in the laboratory polar coordinate are denoted as O, (p, while the polar angle and azimuthal 
angle in the magnetic polar coordinate as 9, tp. Similar to Lee et al. (2010), we use bold type to label the vector or the matrix; while we use 
subscripts XlVlZ and x >. y \ z t to indicate their components in the laboratory frame and the magnetic coordinate, respectively. 

2.2 Geometry, emission region and modeling 

We assume that pulsars have a dipole magnetic field in the magnetic nonrotating frame. Thus dipole field line function can be expressed as 

r = R c sin 2 6, (1) 

where r is the polar radius from pulsar center to an emission point; R c is the maximum radius of the field line, which is the function of a and 
tp s . To obtain the field line constant R c for the last open field line, Zhang et al. (2007) and Lee et al. (2009) derived a cubic equation, i.e., 

A cot 3 M + B cot 2 9m + C cot 6 M + D = 0, (2) 

where A = 4 sin 2 a, B = 5 sin 2a cos -0 S , C = 4 sin 2 a — 6(sin 2 a cos 2 ip s — cos 2 a), D — — sin 2a cos ips, 0m is the polar angle at point 
M. Therefore cot 9m can be determined analytically, and R e is 

D , , s E LC (l + cot 2 6' AJ ) 3 / 2 

V 1 + cot 9m — (cos a cot 9m — sin a cos ip B ) 2 

where i?LC is the radius of light cylinder. Given the values of R e and the last open magnetic field line can be uniquely defined. The polar 
angle 9n of the null charge surface (defined as tt ■ B = 0) is given by 

1 . \/ 8 cot 2 a sec 2 ib B +9 ± cot 2 a sec 2 ips \ 

N = - arccos ( — -5 . (4) 

2 3 cot^ a sec^ ip s + 3 

For a critical field line, which is across the intersection of the null charge surface and light cylinder, there is a relation between On and 
9n (Gangadhara 2004, Wang et al. 2006), i.e. 

On = arccos(cos a cos #n — sin a sin #n cos ip s ) . (5) 
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Then, the critical field line constant, R c ,n, is given by 

i? 0jN (a, ip s ) = -Rlc esc 2 9 N esc 9 N . (6) 

The height of the null charge surface on the last open field line can be derived, i.e. r^(ip B ) = R c (a, ip B ) sin 2 9^. Qiao et al. (2007) gave 
an one-dimension solution to the acceleration electric field, the true 3D solution for the Possion equation with mixed boundary conditions is 
unknown. 

Dyks & Harding (2004) found that, at low altitudes, a distorted angle of the order of (r/R^c) 2 is attributed to the rotation deflection on 
the local direction of the magnetic field, where r is the radial distance of the emission source. In the annular gap model, the radiation region 
is concentrated on the vicinity of the null charge surface of the last open field lines, the emission height is well below the light cylinder radius 
Rlc, and the sweepback effect on open field lines could be ignored. We use the 3D vacuum static dipolar field in this paper. Adopting the 
method of the open volume coordinates (Cheng et al. 2000, Dyks & Harding 2004, Tang et al. 2008, Harding et al. 2008), we calculate the 
polar shape of the annular gap. The inner edge and outer edge of the annular gap region are defined as the footpoints of the critical field lines 
and last open field lines. Here we adopt the conventional wisdom and assume that the 7-ray emissivities I(9 B , ip B ) on each open field line has 
a Gaussian distribution, i.e., 

rta 1 \ t (a i\ r ( L ( s *> Vfe) - L (6 pBak , ip B )) 2 . 

7(0 8 , ip s ) = 7p eak (0p ea k, ipa) exp [ 2o-2 ], (7) 



where L(9 B ,ip B ) = L B \/r 2 + (dr/d9) 2 d9 is the arc length of the emission point on each field line counted from the pulsar center, a is a 
bunch scale of the emission region on each open field line in the annular gap, and Lo(0 pca k, ips) is the arc length of the peak emissivity spot 
P($ P cak, ips) on this open field line. In principle, the P(0 pC ak, ipa) could be located at anywhere on this field line. However, we found that 
our 1-D solution to the acceleration electric potential of each open field lines in the annular gap reaches the maximum near the null charge 
surface. The peak emissivity of the charged particles accelerated by the magnetospheric electric field then should be located near the null 
charge surface, which is also proved by our simulated light curves (See figures in §3). 

The dominated emission region is located near the null charge surface (Qiao et al. 2004). The heights r pca k for the peak emissivity spot 
on each open field line (with the same ip s , different polar angle 9) in the annular gap can be written as 

r P eak(^s) = Akxn(V's) + (1 - A)ktn(O), (8) 

where k is a model parameter for the ratio of the peak emission height with respective to the null charge surface height rt^(ip a ); A is another 
model parameter, describing the deformation of emission location from a circle (Lee et al. 2006); ?*n(0) is the height of the point with 
magnetic azimuthal ip — 0° on the last open field line. Then the peak emission position P on each open field line can be uniquely determined, 
i.e., pC ak = arcsin(\/r poa k/i? Ci f (a, ips)), where R c> { (a, ip s ) is the field line constant of the open field line with ip B . 

The peak emissivity 7 pea k(# P eak, ips) for different open field lines could follow a Gaussian distribution (Cheng et al. 2000, Dyks & 
Rudak 2003), i.e., 

r i ft ni,\-r e ™ r (ds P (ip s ) - 8 cp (ip s )) 2 

J p e a k(fpeak, Ips ) — 10 CXp [ ^ J, (9) 

^ P cak 

where Iq is a scaled emissivity, 9 Bp is used to label a field line in the pulsar annular regions, 9 cp — (0n,^ s + P ,^ a )/2 is the central field 
line among those field lines with ip s . During simulations, we take the width <7 pe ak ~ 0.002 for young pulsars, and ~ 0.008 for millisecond 
pulsars. 

To calculate the light curves for a pulsar, we divide the polar shape of the annular gap into 3 1 rings and obtain the open field line with ip s 
and 9 B . Then we calculate the emissivities projected onto the sky I (9 j, ipj) and direction na (9, , ipj ) of the emission spot (9j , ipj ) on each open 
field line of each ring in the magnetic frame. We also take the aberration effect into account, and use the Lorentz transformation matrix to 
transform the emission direction n.B(9j, ipj) to the direction (<^>j , Cj ) m me l arj f rame (observer frame), where <pj — arctan(n yi y/n I/iX ) is 
the emission spot's rotation phase with respect to the pulsar rotation axis, and (j = arccos(n„, z / ^Jn v ^ 2 + n V:Y 2 + n v ^ 2 ) the viewing angle 
for a distant, nonrotating observer (See details for aberration effect in Lee et al. 2010). We also add the phase shift S(f> rct = r n cos(0 M j — 0j), 
the first order of equation (33) in Gangadhara (2005), caused by the retardation effect to <pj, i.e. cp'i = arctan(7i„,y/n„ lX ) — <50rct, where 
r n = t/Rlc is the emission radius in units of the light cylinder radius and 0,,j is the half opening angle of the emission beam at the emission 
spot {9- s , ipj). 

The "photon sky-map", and the corresponding light curve cut by a line of sight with a viewing angle £ are therefore finally derived. The 
application of the annular gap to young and millisecond pulsars for high energy pulse profiles are presented in Section 3. 



3 SIMULATED LIGHT CURVES FOR YOUNG AND MILLISECOND PULSARS 

Young 7-ray pulsars are very energetic, due to their large spin-down luminosity. The flux of gamma ray emission detected by EGRET are 
very high, e.g., 10~ 10 — 10~ 8 ergs -1 s~ 2 (Thompson et al. 1999). We use the annular gap model to simulate the light curves for the three 
brightest pulsars: the Crab, the Vela and the Geminga, respectively. The simulated light curves can reproduce the main observation features 
detected by Fermi (Abdo et al. 2010b). Most of the distinct light curves of millisecond pulsars detected by Fermi (Abdo et al. 2009a) can also 
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Table 1. The parameters for the simulated light curves for young and millisecond pulsars. 
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Figure 2. The 7-ray photon sky-map (observer angle f vs. rotation phase <f>) and simulated light curves for the Crab pulsar (left panels), the Vela pulsar 
(middle panels) and the Geminga pulsar (right panels) in the framework of single-pole anuular gap model, in comparison with the observations (thin solid 
lines, taken from Figure A-8, A-16 and A-12 of Abdo et al. 2010b). The simulated light curves are smoothed to 64 bins. The phase of magnetic pole in the 
photon sky-maps is set to be -0.15 to show the light curves in the middle. The observed light curves (thin solid lines) are compared with the simulated ones 
(thick solid lines). Parameters are given in Table 1. 



be explained in the annular gap model, especially for pulse profiles of a single peak or small separation (~ 0.2 in phase). We have simulated 
three light curves of millisecond pulsars: PSR J0030+0451, PSR J0218+4232 and PSR J0437-4715, respectively. The relevant results will be 
briefly described below. The related parameters from simulations are given in Table 1. 



3.1 The Crab pulsar 

The Crab pulsar is famous for the phase-aligned pulse profiles of multi-wave band emission and beautiful X-ray torus. The observed 7-ray 
(> 100 MeV) features of the Crab pulsar (Abdo et al. 2010a and Abdo et al. 2010b) are two peaks with a separation of ~ 0.4 rotation phase 
and a hard bridge. The observed light curves at different bands are quite similar. We simulated multi-band light curves for the Crab pulsar, but 
presented only one band in Figure 2 (left panels). It shows most of the observed features, especially for the phase-aligned pulse profiles with 
nearly similar peak ratio (P2/P1). An important result is that the 7-ray emission with higher photon energy comes from higher altitudes in the 
magnetosphere. This is probably caused by strong 7-B absorption effect for emission from lower altitudes (Lee et al. 2010). In addition, the 
parameters of a and £ in our model are based on the adopted value (Wang 2003, Harding et al. 2008) and the results of X-ray torus simulation 
(Ng & Romani 2008). Other best k and a values indicate that the 7-ray photons are emitted from a wide range of emission altitudes from the 
neutron star surface to about half of the light cylinder radius. 



3.2 The Vela pulsar 

Similar to the Crab pulsar, the Vela pulsar is the brightest object in the 7-ray sky. Figure 2 (middle panels) shows the simulated light curve 
of > 0.1 GeV band for the Vela pulsar. The main features, two sharp peaks with a separation of ~ 0.42 and the peak ratio observed by Fermi 
and EGRET (Thompson 2001, Abdo et al. 2010b), are approximately reproduced. Similarly, we choose the parameters £ according to the 
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Figure 3. The same as Figure 2 but for the millisecond pulsar PSR J0030+0451 (left panels), PSR J02 18+4232 (middle panels) and PSR J0437-4715 (right 
panels) in comparison with the observations (thin solid lines, taken from Figure A-2, A-4 and A-7 of Abdo et al. 2010b, respectively). The simulated light 
curves (thick solid lines) are smoothed to 32 bins. 



result of X-ray torus simulation (Ng & Romani 2008). The inclination angle a = 30° which is thought to be an intermediate inclined rotator 
for the Vela gives the "best" simulated results. The best parameters k and a indicate that the 7-ray photons mainly come from high altitudes. 
Here we only used the open field lines in the annular gap region for the simulations. To get better results, especially for the third peak in the 
bridge, both the annular gap and core gap regions (Qiao et al, 2007) probably should be used to simulate the observed light curves for the 
Vela pulsar. 



3.3 The Genminga pulsar 

The Geminga pulsar is another bright 7-ray pulsar. Figure 2 (right panels) shows the simulated light curve for Geminga. The peak separation 
and peak ratio (Abdo et al. 2010b) are reproduced. The best k and a values indicate that the emission region could be mainly above the 
height of the null charge surface. 



3.4 PSRJ0030+0451 

PSR J0030+0451 is a recently discovered solitary milisecond pulsar. The pulsed 7-ray emission was detected by Fermi (Abdo et al. 2009b). 
The simulated light curve (> 0.1 GeV) for PSR J0030+0451 is shown in Figure 3 (left panels), which is similar to the observed features. The 
parameters a is chosen for small magnetic inclination angles (Zhang et al. 1998, Tauris & Manchester 1998, Young et al. 2009). Furthermore, 
the k and a values indicate a small emission region at low altitudes. 



3.5 PSRJ0218+4232 

The millisecond pulsar PSR J0218+4232 has a complex observed pulse profile (Figure A-4 of Abdo et al. 2010b). Figure 3 (middle panels) 
shows the simulated light curve (> 0. 1 GeV) for PSR J02 18+4232. Our result roughly reproduces the observed peaks. The model parameters, 
similar to those of PSR J0030+0451, also favor a small emission region at low altitudes. 



3.6 PSRJ0437-4715 

PSR J0437-4715 is the nearest millisecond pulsar with a good radio timing. It has a single narrow (about 0.2 phase) pulse peak at 7-ray 
band (Abdo et al. 2010b). Figure 3 (right panels) shows the simulated light curve for PSR J0437-4715. Our simulated result reproduces the 
features with a small inclination angle a and a large viewing angle (. Again, the model parameters are similar to those of PSR J0030-0451 
and also favor a small emission region at low altitudes. Note that the 7-ray emission beams (photon sky-maps) of millisecond pulsars have a 
hollow cone in shape. 
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4 DISCUSSIONS AND CONCLUSIONS 

As an approximation, a static dipole field is used to model the pulsar magnetosphere. The pulsar magnetic field can be approximated by a 
static magnetic dipole configuration if the radial distance is not so far away from the pulsar surface (Muslimov & Harding 2005). The critical 
magnetic field lines used to define the annular gap are different in the magnetosphere models such as the retarded vacuum (Cheng et al. 2000) 
and force-free models (Spitkovsky 2006), since the positions of the null surface are different. However, in our model, the 7-ray emission 
regions are concentrated in the middle field lines of the annular gap, and the peak emission comes from the vicinity of the null charge surface. 
The high energy emission from the field lines near the upper boundary (critical field lines) and the ones near the lower boundary (last open 
field lines) give little contribution to the observed light curves. Our model is different from caustic models, e.g., the outer gap model, the slot 
gap model and the separatrix layer model (Bai & Spitkovsky 2009b), which all assumed a uniform emissivity along a field line. Our simulated 
light curves are mainly dependent on the non-uniform emissivities in a deformed radiation beam (see A in equation 9) of the annular gap 
region. A is described in detail in Lee et al, (2006), the large value leads to a more deformed radiation beam from the circle one. Owing to the 
intermediate height of 7-ray emission region, our annular gap model is weakly dependent on the magnetic field configuration either static or 
retarded dipole field. 

The annular gap has a sufficient thickness of trans-field lines and high altitude acceleration regions. This leads to a fan-beam 7-ray 
emission and sufficient photon luminosity, which is suitable to interpret the observed light curves and the broad-band emission. The annular 
gap model holds the advantages of the slot gap and outer gap models, and works for pulsars with short spin periods. We use the annular gap 
model to simulate the light curves for young and millisecond pulsars with three assumptions: (i) the emissivities on a single field line between 
the critical field line and the last open field line follows a Gaussian distribution; (ii) the peak emission spot of a single field line is located at 
the vicinity of the null charge surface; (iii) the peak emissivities of a group of field lines with the same magnetic azimuthal (i.e., in the same 
plane) between the critical field line and the last open field line follows another Gaussian distribution. The assumption (ii) is consistent with 
our 1-D solution to the acceleration electric potential drop in the annular gap. The other two assumptions are mainly based on the magnetic 
pair absorption and the 3D global parallel electric field for the acceleration of relativistic charged particles in the magnetosphere, which will 
be investigated in future. In our calculations, it is shown that the emission region extends from the neutron star surface to about the half of 
the light cylinder radius, and our annular gap model is an intermediate emission height model. We also find the following conclusions from 
the modeling. 

(1) The simulated light curves can reproduce most of observed features for both young and millisecond pulsars. The 7-ray emission 
with higher photon energy comes from higher altitudes in the magnetosphere. 

(2) The 7-ray beams for both young and millisecond pulsars are hollow cones in shape. 

(3) The 7-ray emission light curves (pulse profiles) are determined by the inclination angle of magnetic dipole field and the observer's 
viewing geometry. The magnetic inclination angles and viewing angles of millisecond pulsars can not well constrained by any methods at 
present. This leads to difficulties for precise reproduction of light curves. 

(4) The radiation regions of young pulsars are larger. In terms of an individual field line, this large region covers from the pulsar surface 
to about the half of the light cylinder radius. The peak emission comes near the null charge surface. 

(5) The radiation regions of millisecond pulsars are small, from the pulsar surface to about the one third of the light cylinder radius. The 
peak emission comes from a region below the null charge surface. 

(6) Our model favors small inclination angles (a < 35°) for the millisecond pulsars, and larger inclination angles (a ~ 30° — 70°) for 
the young pulsars. This is somewhat compatible with the alignment of the spin and magnetic axes over about 1 Myr from the analysis of the 
new pulsewidth data (Young et al. 2009). 

(7) Our results also show that the solid angle of gamma-ray beams are much less than 4-7T, especially for MSPs. This will reduce the 
7-ray emission conversion efficiency r\ = L^/Eyot, and solve the puzzle for 77 > 1 in Abdo et al. (2010b). 

As shown in Qiao et al. (2004; 2007), the annular gap can have sufficient electric potential drop to produce pairs that can generate radio 
emission. Radio emission can be generated even at a higher region either for the annular gap or core gap. If radio emission comes from a 
lower region of either inner annular gap or the core gap (Qiao et al. 2004), the radio peak should appear between the two 7-ray peaks. On the 
other hand, it is possible that the radio radiation comes from a higher region or the opposite magnetic pole, producing a leading or trailing 
radio peak. 
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